Potential and field singularity at a surface point charge 
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The behavior of the magnetic potential near a point charge (fluxon) located at a curved regular 
boundary surface is shown to be essentially different from that of a volume point charge. In addition 
to the usual inverse distance singularity, two singular terms are generally present. The first of them, 
a logarithmic one, is axially symmetric with respect to the boundary normal at the charge location, 
and proportional to the sum of the two principal curvatures of the boundary surface at this point, 
that is, to the local mean curvature. The second term is asymmetric and proportional to the 
difference of the two principal curvatures in question; it is also bounded at the charge location. 
Both terms vanish, apparently, if the charge is at a planar point of the boundary, and only in this 
case. The field in the charge vicinity behaves accordingly, featuring generally two singular terms 
proportional to the inverse distance, in addition to the main inverse distance squared singularity. 
This result is significant, in particular, for studying the interaction of magnetic vortices in type II 
superconductors. 
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Magnetic vortex lines are formed in type II 
superconductors^. When crossing the superconductor 
boundary, they create strongly localized surface sources 
of magnetic field (fluxons), which may play an impor- 
tant role in various physical situations. For instance, 
two space tests of Einstein's General Relativity, Gravity 
Probe B2* 3 . and STEP (Space Test of the Equivalence 
Principle)—, are based on low temperature technology 
with type II superconductors, and their setup is signifi- 
cantly affected by fluxons. 

The size of a surface magnetic spot is about the micro- 
scopic London lengthi, i. e., it is typically much smaller 
than characteristic macroscopic sizes involved. Thus the 
point charge approximation appears naturally and proves 
to be sufficient for many applications. Within this ap- 
proximation, the magnetic potential, ip = i/j(R), satisfies 
the Neumann boundary value problem 



If D is bounded, then each vortex line starts and ends 
at the boundary, the number of charges is even, and the 



total charge vanishes, ^ 



N 



1 Vj *o 



0, which condition is 



dtp 
On 



Avb = 0, Refl, (1) 

N 

= v i $ o S S (R - Rj)) R, Rj e S . (2) 



Here the domain D is the empty space, surface S is 
the superconductor boundary, $0 = h/2e is the mag- 
netic flux quantum in SI units, and the magnetic field is 
B = —Vip. Moreover, <$s(R — Rj) denotes the surface 
delta- function at the position, Rj, of a vortex, and Vj is 
either plus or minus one, depending on whether the field 
line enters the domain D (vj = +1), or exits it (yj = —1). 
We assume that the boundary S is smooth enough (at 
least C 3 ) near every charge. Outside the charge vicini- 
ties it may have any singularities compatible with the 
finite local energy condition, meaning (V"0) 2 is locally 
integrable. 



the solvability criterion of the problem, Eqs. ^0). If the 
domain D is infinite, some field lines may end at infinity, 
and this condition may not hold; in any case, we do not 
use it in the following analysis, which is entirely local. 

An immediate question regarding the above boundary 
value problem is how does its solution behave near a sur- 
face charge? For a curved boundary, an answer based on 
the similarity with the volume point charge turns out in- 
correct. This is seen from the simplest example, a spher- 
ical domain. A closed-form exact solution to Eqs. fflEJ) 
in the exterior of a sphere was obtained in Ref. It 
shows that a new logarithmic singular term, inversely 
proportional to the radius of the sphere, is added to the 
main inverse distance singularity in the expansion of the 
potential near the charge. So, what happens with the 
singularity for a generally curved smooth surface? 

Our search for the answer to this natural and, in 
fact, classical question covered books and papers in both 
mathematical physics and in the field of vortices in super- 
conductors, as well as communications with colleagues in 
both fields. We also talked with high energy theorists ex- 
pecting to find perhaps some relevant results in view of 
the discussions of the magnetic monopole. However, no 
ready answer was found, which might be not so surpris- 
ing. Indeed, the Neumann boundary value problem with 
surface charges is not relevant to the design of electro- 
static systems. On the other hand, its magnetostatic im- 
plementation became available only with the widespread 
technical use of superconductors in the recent years. Last 
but not least, the answer proves to be not that simple. 
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superconducting bulk 



F(x,y) 




empty space, D 



FIG. 1: Definition of coordinate systems near a charge. 



In this letter we fill the gap by deriving a complete sin- 
gular part of the expansion of the solution to Eqs. (^Q [2I> 
near a charge at an arbitrary curved smooth boundary. 
As compared to the case of a sphere, one more singular 
term, proportional to the difference of the two principal 
curvatures, appears in the general case. 

We are interested in the behavior of the potential near 
a single surface charge at some R 3 . For brevity, we thus 
drop the charge index in the following calculation. We 
put the origin of a Cartesian coordinate system at Rj, 



so that r 



R — Ro 



We point the z axis along the 



outward normal to the surface S (that is, into the su- 
perconducting bulk), choosing the x and y axes in the 
tangent plane, so that the unit vectors {x, y, z} form a 
right orthogonal triplet. Along with Cartesian {a;, y, z}, 
we will use the corresponding spherical, {r, 9, </>}, and 
cylindrical, {p, (p, z}, coordinate systems (see Fig.0. 

The shape of the smooth boundary surface in the 
vicinity of the charge can be described by the equa- 
tion z — F(x,y). The Taylor expansion of the function 
F(x, y) around x = y — apparently has no terms linear 
in x or y, since z is oriented along the normal. More- 
over, by an appropriate rotation of the coordinate axes 
x, y in the tangent plane, we can ensure that the second 
cross-derivative of F vanishes at the origin, hence the 
expansion acquires the form 

z = F(x,y) = — x 2 + — y 2 + 0(p 3 ) = f(x,y) + 0(p 3 ) , 

(3) 



where 



fc<*> = 



d 2 F 



dx 2 



k (y) = 



d 2 F 



x—y—0 



dy 2 



(4) 



x—y—0 



are the two principal curvatures of the boundary surface 
at the charge location. 



Since, near r = 0, 
5 s (r) = 5(x)5(y)/J, 

d/dn = n- V= (l/J)(d/dz-F x d/dx-F y d/dy) ; 
J = Jl + FZ + F*, 

the boundary condition, Eq. (0), in the vicinity of the 
charge can be written in terms of variables x, y, z &si&: 



dtp 



is<5>o$(x)6{y) 



z=F{x,y) 



F, 



dtp 
dx 



F 



dtp 
dy 



(5) 



z=F(x,y) 



The partial derivatives of the function F(x, y) near the 
origin are given, to the order we are interested in, by 



F x — k^x + 0(p 2 ), F y = k^y + 0(p 2 ) 



(6) 



Once again, we only care about the vicinity of the 
charge where z = F(x, y) is small, so we can use pertur- 
bation of the boundary to move the boundary condition, 
Eq. JSJ, to the plane z = 0. This is done by means of 
the following Taylor expansion of an arbitrary function 
w = w(x, y, z): 



it' 



z=F(x,y) 



+F 



dw 

dz 



F 2 d 2 w 



z=0 



2 dz 2 



2=0 



Applying this to the derivatives of ip m Eq. (J5J we write 
it, to the proper order, in the form: 



dtp 
dz 



= v<f> 8(x)5(y) 



dip 
dx 



p d± p d 2 tp 
y dy dz 2 



(7) 



The final step of this derivation is to expand iJj in a 
series of successively smaller (that is, less singular at the 
origin) functions tp , 



ijj = ?a (0) + ip {1) + tp {2) + 



(8) 



Introducing this expansion in the boundary condition 
Eq. 10 and then matching the terms of the same or- 
der, we end up with the following sequence of boundary 
conditions for tp^\ i = 0, 1, 2, . . . , at z = 0: 



dz 



v$ 5(x)6(y) ; 



dz 



z=0 



z=0 



c^ (0) , ,.M..&i> l0) 



-f{x,y) 



dx 

52^(0) 
dz 2 



k (y) 



dy 



(9) 



(10) 



2 = 



and so on. Here we have dropped higher order terms in 
the right-hand sides by replacing F and F x , F y with their 
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main term expressions from Eqs. © and ©, respectively. 
Of course, all functions rpw are subject to the Laplace 
equation, Eq. @J. 

Thus, locally we have successfully replaced the bound- 
ary value problem of Eqs. Q 0) in the domain D by a 
sequence of problems for functions ip^\ i — 0, 1, 2, ... , 
harmonic in the half-space z < and satisfying the above 
boundary conditions, Eqs. ©, l(TU|> . etc. We now need to 
solve these problems for the half-space one by one, until 
the normal derivative of the solution becomes finite at 
the boundary. 

The zero order solution ip^ obeying the boundary con- 
dition of Eq. 10 is, of course, 



(0) 



i/$ 1 

27r r 



(11) 



It allows one to immediately calculate the r. h. s. 
Eq Indeed, 



of 



f 



#yo) 



dz 2 



^1 



v$of ( 1 
j/$ k^x 2 + k^y 2 

~27T V 



where the second term in the middle expression turns to 
zero at z — 0, contributing no 5-like singularities, due to 
the presence of the factor / = 0(p 2 ). Taking also into 
account that 9r _1 jdx — —x/r 3 , dr^ 1 jdy = —y/r 3 , we 
find the boundary condition for ipW in its final explicit 
form: 



dz 



z/$ k^x 2 + k^y 2 



2tt 

_ fc(w) 
P 



2p 3 

+ fete) 



■ cos 2<p 



(12) 



The two terms on the utmost right here have essentially 
different singularities at the origin. For this reason, we 
treat them separately by splitting the problem in two in 
the following way: 



(1) = 



dtp. 



(i) 



dz 



dz 



8^ p 

u$ k( x) - fc(») 



8tt 



P 



cos 2(f> 



(13) 
(14) 

(15) 



The Neumann problem for ip^ in the half-space does 
not have solutions bounded at infinity, as one would ex- 
pect in our investigation (we are actually looking for 
terms growing away from the charge, because a weaker 
singularity next to the inverse distance is most probably 
some logarithm tending to infinity at both the charge and 
the infinite distance from it). For this reason, no solution 



can be found by means of standard techniques. However, 
a harmonic and regular in the half-space z < function 



(i) 



K + \n[{r - z)/d] 
if + [ln(r/d)+m(l-cos6»)] 



K+ = V<frn 



k ix) ± k {v) 



(16) 
(17) 



where d > is an arbitrary constant of the dimension of 
length, provides the needed solution. Indeed, it satisfies 
the boundary condition, Eq. I|14|). in view of 



<91n(r — z) z/r — 1 

dz r — z 



-l/r->-l/p, 



-0 . 



The solution given by Eq. (|16|l is unique in the class of 
functions with the logarithmic growth at infinity, namely, 
those with the asymptotics 

ipW = K + ln(r/d) + K + \n(l-cos6) + o(l), 
d^ 



di 



K+/r + 0(l/r 2 ), 



oo 



Contrary to the previous one, the Neumann problem 



for Ipr^, 

A^ 1} =0, z < 0, 



(1) 



dz 



z=0 



cos 2(f> 

P 



,( 18 ) 

has a unique, up to a constant, solution bounded at infin- 
ity [namely, a solution that obeys somewhat unusual con- 
ditions ip^ = O(l), dip^/dr = o(l/r 2 ), r — » oo]. 
The solution is obtained by the standard separation 
of variables in cylindrical coordinates using the Hankcl 
transform, and it reads: 



(i) 



K- cos 20 
K- cos 2<p 



00 dX 
J 2 (Ap)exp(-A|z|) — 
ii A 

2 

P 



~2~ 



y 



(r - zf 



(19) 



The value of the integral is found in Ref. |g, 4.14.(5), and 
the constant K- is defined in Eq. p7JI. Interestingly, this 
solution in spherical coordinates does not depend on the 
radius, being a function of the angles only [singular on 

the positive semi-axis z > 0, same as ip^ in Eq. O]: 



4 1] = ~ 



K_ sin 2 9 cos 2<p 



dip- 



in 



2 (l-cos(9) : 



dr 







It is now straightforward to see that the Neumann 
boundary data for all higher order corrections to the po- 
tential, starting with ip^ 2 \ are finite at the origin (and 
dropping fast enough at infinity); accordingly, the solu- 
tions of the corresponding problems bounded at infinity 
are unique up to an additive constant. It also means that 
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all the terms in the expansion, Eq. (JSJ , of the potential, 
whose normal derivative are singular at the location of a 
surface charge, are given by the solutions already found. 
Hence, combining the expressions from Eqs. (|16H . 
and (|19|l . we find the desired formula for the magneto- 
static potential near a surface charge (r — » 0): 



Eq. l(2Tj|l becomes 



v> (0) + 4 1] 

2ir 



(i) 



1 fcW + fe(v) , r - 

- + * ln — T 

r 4 d 



fc(x) _ fc (y) x 2 _ y 



2 1 



(r — z) 2 



(nonsingular 
terms) 



(20) 



It is easy to rewrite this in our general notations from 
Eqs. HI I2J) by replacing |r| with |R — Rj|, x with X-Xj, 
etc. Instead, we give the expression of the singular part 
of the magnetic field near the charge. It can be written 
in the form: 



B 



2ir 



4r 



fcO<0 _ fe(») 



4r (l-cos(9) 2 
-(nonsingular terms) 



sin 6* 

1 — cos t 

cos2cj)9 + sin 20 1 



(21) 



Here are a few concluding remarks regarding the ob- 
tained result. 

First, notice that the leading order contribution to the 
potential, Eq. (|llfl . is twice that of the point charge lo- 
cated in a volume away from its boundaries. This is 
clearly explained by the fact that the field lines and the 
flux from the surface charge emanate only into the half- 
space, versus the full space for the volume charge. 

The two singular corrections to the usual inverse dis- 
tance singularity of the potential, Eq. I|2(J|) . are very dif- 
ferent. The first one is logarithmic, symmetric about the 
direction of the normal to the boundary at the charge 
location, and proportional to the sum of two principal 
surface curvatures there, i. e., to the mean boundary cur- 
vature. Thus, it vanishes if the charge sits at a symmet- 
ric saddle point of the boundary. The second additional 
singularity is asymmetric, proportional to the difference 
of the principal curvatures, and vanishes thus when the 
latter are equal, i. e., when the charge is at a spherical 
point of the boundary. This second term is bounded at 
the charge location [giving unbounded field components, 
see Eq. I|21|)]. but is not uniquely defined there, with the 
limiting values depending on the direction along which 
the limit is taken. Note that both corrections vanish si- 
multaneously if and only if the charge is at the planar 
point of the boundary. 

In a particular case when the domain D is the exterior 
of a sphere of the radius a, one has k^ — = I /a. If 
there is just one surface charge, N = j = 1 and v = 1 
(so that the incoming vortex line ends at infinity), the 



VJ= 2^ 



1 



1 IR-RJ -n- (R-Ri) 
+ 2^ ln d 



|R-Ri 

+ (nonconstant nonsingular terms) , (22) 

in complete agreement with the exact solution obtained 
in Ref. |5| with d = 2a. 

Finally, the obtained singular expansion of the poten- 
tial can be used in the derivation of the force acting on 
a charge in a fashion similar to the one developed in 
the case of volume point charges^, i. e., by means of the 
geometrical regularization of energy and, henceforth, the 
force, as the energy gradient in the charge location. How- 
ever, in a striking contrast with the volume case, the force 
here is found to depend on the gradient of the curvature 
at the charge location. Namely, due to the first additional 
singular term in the potential, Eq. I|2U|) . there appears a 
tangential force on the charge which tries to move it to- 
wards the point of the stationary mean curvature of the 
boundary, and which diverges in the regularization limit. 
If confirmed, this divergence would mean that cither the 
approximation of the point surface charges does not com- 
pletely describe real microscopic, but finite size fluxons, 
or, strangely enough, that the fluxons cannot reside at 
arbitrary points of a curved boundary, or perhaps even 
something else. 

A detailed study of fluxon interactions will be carried 
out in a separate publication. However, it is clear that it 
will necessarily use the results of this paper, in view of 
the relation 



Js 



(23) 



where ^(r, A*) * s the potential created by the surface 
charge density £ &i an d VK r ) i s the potential 

from Eq. I|20|l . For small, yet finite size fluxons the diver- 
gent asymptotics derived above will have an explicit short 
scale cutoff defined by the spatial extent of the density 
(presumably, the London length). However, the detailed 
analysis will require a deeper insight in the real struc- 
ture of magnetic vortex lines near a boundary. Without 
such an analysis one cannot, in fact, speculate about the 
strength and importance of these surface interactions; we 
will thus limit ourselves to just a few short comments. 

First, one compares, naturally, the surface force com- 
ing from the logarithmic term in the field potential to 
the strength of the random pinning force that defines 
the fluxon's position^. The latter depends on the flux 
tube length and the former does not. So, allowing for a 
physical regularization of the mathematically divergent 
surface effects, one will in any case come up with some 
characteristic length, L, below which the surface force 
will dominate. The description of the vortex line dynam- 
ics that does not account for surface effects at distances 
from the surface smaller than L is necessarily incomplete. 

Second, forces between two vortices in a superconduct- 
ing bulk are exponentially small if the vortex line sepa- 
ration is larger than the London length (precisely the 
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regime we are discussing) . These forces can be neglected. 
Thus, the surface effects we have found will be the leading 
interaction terms. Such effects are significant and trans- 
late into an experimentally relevant magnetic "friction" 
between superconducting bodies 9 . 
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